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The phase diagram near the magnetic field-induced quantum critical point of dimerized spin systems can 
be extremely rich when the interactions are frustrated. Complex phases that emerge near this critical point 
are derived under control by mapping the spin system into a dilute gas of bosons. Geometric frustration of 
the dimer lattice manifests in a rather large degeneracy of wave-vectors Q^^ that minimize the single-boson 
dispersion. The bosons can then condense into a single or a multi-Q state. The magnetic counterparts of these 
condensates range from a simple spiral (single-Q condensate) to magnetic vortex crystals (6-Q condensates). 
Inspired by the *S = 1 dimer antiferromagnet Ba3Mn208, we calculate the eff'ective boson-boson interactions in 
the dilute limit and demonstrate that a plethora of exotic condensates emerge above the critical field H = Hd 
and that a six-fold symmetric exchange anisotropy selects a vortex crystal or a double-Q fan state. 

PACS numbers: 



The emergence of topological spin textures in solids trig- 
gered an enormous interest because of their relevance for spin- 
electronic technology. Outstanding examples are the magnetic 
skyrmion crystals that have been recently observed in spe- 
cific alloys with B20 structure.^ ^ The helimagnetic ground 
state of these materials results from competition between fer- 
romagnetic exchange and Dzyaloshinskii-Moriya (DM) inter- 
actions. The skyrmion lattice emerges in a narrow temperature 
and magnetic field window just below the magnetic ordering 
temperature.^ Last year, Seki et ah reported the first discov- 
ery of a spontaneous magnetic skyrmion crystal in the insu- 
lating chiral magnet Cu20Se03^ and opened a new promising 
route to electric control of magnetism. Helical magnetic or- 
derings, such as skyrmion lattices, induce another texture of 
electric dipole moments enabling a magneto-electric response 
that can be exploited to manipulate the magnetic state with 
electric field gradients.'^ This procedure is energetically more 
eflftcient than current-driven approaches required for metals. 

After this sequence of discoveries, it is natural to ask if 
topological spin textures can emerge under more general con- 
ditions. In particular, we would like to know if these textures 
can arise from frustration between exchange interactions in- 
stead of competition between the exchange and DM coupling. 
The advantage of the first mechanism is that the ratio between 
exchange constants can be changed with pressure, while the 
ratio between exchange and DM couplings is much less tun- 
able. Thus, the first mechanism enables pressure-induced tun- 
ing of the spin texture structural parameters. To explore this 
exciting possibility, we focus on the field-induced quantum 
critical point (QCP) of a frustrated spin dimer Mott insulator 
that is a quantum paramagnet at if = 0: each dimer is pre- 
dominantly in its singlet state. The first excitation is a single- 
dimer triplet state that propagates with well-defined momen- 
tum (triplon). The QCP signals the onset of a Bose-Einstein 
condensate (BEC) of triplons at the field, H = Hd, that 
closes the single-triplon energy gap.^ While basic concepts 
of this transition were established through the study of ef- 
fectively unfrustrated magnets like TlCuCls,*^ strongly frus- 
trated materials like SrCu2(B03)2p^ BaCuSi206™^ and 
B as Mn20gf^^^ exhibit more exotic behavior. 

Ba3Mn208 consists of weakly-coupled dimers of 5" = 1 



Mn^^ (3J^) ions that form triangular lattice bilayers stacked 
along the c axis in a period- 3 structure (Fig. [T]). Besides 
the intra-bilayer exchange, which would alone favor simple 
three- sublattice commensurate spin orderings,^^ inelastic neu- 
tron scattering (INS) measurements^^ revealed several frus- 
trated inter-bilayer exchange constants that shift the single- 
triplon dispersion minima to six-fold degenerate incommen- 
surate wave- vectors ±Q„ (1 < ^ < 3) [Fig. [TJd)]. While a 
simple single-Q state is stable over a big field interval between 
Hci and Hd (Hd ~ 9 T, Hc2 ~ 26 T), the observation of a dif- 
ferent ordered phase right above Hd suggests the possibility 
of an unconventional multi-Q condensate. 

The distribution of the triplon condensate among the six- 
fold single-triplon dispersion minima is determined by the ef- 
fective triplon-triplon interaction. Here we consider a hard- 
core boson model derived from the spin Hamiltonian of 
Ba3Mn208 for H || c and H < Hd. We use the measured 
single-triplon dispersion, and vary the triplon-triplon inter- 
action parameters that are more difl&cult to measure. By ex- 
panding in the small lattice gas parameter near the BEC-QCP 
we predict a rich phase diagram including two different crys- 
tals of magnetic vortices among other interesting phases. 




FIG. 1: (a) Magnetic (Mn^+) lattice of Ba3Mn208. (b) Nearest- 
neighbor (NN) and (c) next-NN inter-bilayer exchange couplings, (d) 
Single-triplon dispersion. Filled (open) triangles indicate the incom- 
mensurate wave- vectors of the lowest-energy modes. 
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FIG. 2: (a) Effective vertex from ladder diagrams. The wavy line 
(filled square) represents the bare (eff'ective) potential, (b)-(f) Eff'ec- 
tive interactions in the Ginzburg-Landau theory. 



FRUSTRATED HARD-CORE BOSON MODEL 

The S - \ spin Hamiltonian of Ba3Mn208 is defined on 
the lattice depicted in Fig.[lJ^^ The intra-dimer antiferromag- 
netic exchange, /q ~ 19.05(4) K/^ is much bigger than any 
inter-dimer exchange. The single-dimer low-energy sector is 
conveniently described with bosonic bond operators P^l^ that 
create singlet, >s'J|0)/, and triplet states, ^J^|0)/ {cr - 1,0, T), in 

each dimer / and satisfy the constraint s^.s.-\-Y^^t^.t. =1. 
Spin-2 states can be ignored near H - Hd because they 
do not modify the triplon dispersion to lowest order in the 
inter-dimer exchange. At the mean field level, the paramag- 
netic state for H < Hd i^ described as a singlet condensate: 
s^. = s. = \ The resulting Hamiltonian has a quadratic 
contribution in the triplet operators that is diagonalized in mo- 
mentum space, = N'^I^Yji exp(/k • r^) t]^ (N is the num- 
ber of dimers), by a standard Bogoliubov transformation;^^ 



= ^l^l-yl - g/^B^^ ^k = /o - ^ + Tk. (1) 



7^2 includes only the lowest energy = I single-triplon 
branch: is the 5"^ = 1 Bogoliubov quasi-particle (dressed 
triplon) that becomes soft at if = Hd and <jJ\^ is the 
single-triplon dispersion, = (8/3)[(/2 - Jd)?^^^ + ('^i,p - 
^i,x)7k^ + {J4,p - Ax)7k^^] is the bare triplon dispersion 
where y^^^ = cos ki -h cos ^2 + cos (^i - ^2), 7^ ^ = cos -h 
cos (^3 - ^1) + cos (^3 - ^2), and y^^"* = cos (-^1 ^2 + ^3) + 
cos (^1 - ^2 + ^3) + cos (^1 -h ^2 - ^3)- The intra-layer triplon 
hopping amplitude is proportional to J2 - J3 ~ 1.318(8) K, 
where J2 and J3 are the nearest-neighbor (NN) intra-bilayer 
exchange couplings along parallel and cross paths, respec- 
tively. We include both parallel- and cross-path couplings in 
the NN (/i,p and /i,x) and next-NN inter-bilayer exchanges 
(/4,p and /4,x) [Figs. [TJb) and [TJc)]. INS results indicate 



Ji,p - Ju ~ -0.685(10) K and /4,p - Ax ~ -0.215(10) K.CHI 
The inter-layer triplon hopping amplitudes are proportional 
to these diflTerences and they are responsible for the in-plane 
shifts of the oji^ minima from the two commensurate three- 
sublattice ordering wave-vectors (2D limit) to the six wave- 
vectors shown in Fig. [TJd): Qi = (2a, a, a), Q2 = (-a,a,0), 
Q3 = (-a, -2a, -a), -Qi, -Q2, and -Q3 with 



a = arccos 



1 Ji- J?> + Ji,p - Jl,x 



2 J2 - J3 -\- 2 (/4^p - /4^x) 



yx 0.924(5). 



(2) 



The boson-boson interaction terms are 



k,k',q 



Vq = 2 (y2ff + ViTk ^ + ^4rk ^) , V2 = (J2 + J3) /2, 

Vi = (/i,p + /i,x) /2, V4 = (/4,p + Ax) /2. (3) 

Here we used that t^i in the low-density limit. The local 
constraint sjs. -h 2o- ^Jo-h a ~ imposed by including an on- 
site hard-core potential U ^ 00}^ Vi, V2, and V4 are the oflf- 
site interactions determined by the exchange constants. In the 
long wave-length limit, k ^ ±Q„, the eflTective triplon-triplon 
interaction is obtained by adding the ladder diagrams shown 
in Fig. [2]^ a). ^^ ^^ The resulting interaction vertex rq(k, k') for 
incoming particles with momenta k and k' and the momentum 
transfer q is asymptotically exact in the dilute limit H ^ Hd . 
(The computation of rq(k, k') is described in Methods [A[) 

GINZBURG-LANDAU THEORY 

If N~^^^{t+qJ = -yJP+Q exp(/0+Qj is the condensate wave 
function of the lowest-energy mode ±Q„ (1 < ^ < 3), the 
triplon density isp = YjniPq +P-q ) 1) and the Ginzburg- 
Landau (GL) expansion of the ground state energy density is 



n n 
n n<m 

+ r4 J]K„P-Q,„ +P-q„PqJ 

n<m 
n<m 

where yu = gfi^(H-Hd) and 0„ = +0-q„- The coeflftcients 
are obtained by replacing the bare interactions in *K4 by the 
eflTective interaction vertices: 

Fi = Fo (Q„ Qn) , F2 = Fo (Q„, -QJ + F_2q„ (Q., -Q.) , 
F3 = Fo (Q„, Q;^) + Fq„^_q„ (Q„, Q^) (n ^ m), 
F4 = Fo (Q„, -Qm) + F-Q„^-Q„ (Q„, -Qm) (n ^ m), 

= Fq„^-q„ (Q„, -Qn) + F_Q„^_Q„ (Q„, -QJ (n ^ m). (5) 
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FIG. 3: (a)-(f) The ground state phase diagram at the BEC-QCP as a function of the bare interaction potentials Vi, V2, and V4. The dashed 
hne indicates a boundary above (below) which < (> 0), while the double line shows a phase boundary between the state No. 2 (below) 
and No. 8 (above) when weak anisotropic perturbations are taken into account (see text), (g) Schematic pictures of the condensate distributions 
among +Q„ (I < n <3). The arrows representing 0„ = 0q„ + 0_q,^ are shown only for those states with particular conditions for 



Most terms in ^ are density-density interactions [Figs.[2];b)- 
[2je)]. Only the T5 term of the ofF-diagonal form 
Q ^-Q ^-Qy. ^ sensitive to the relative phases 

o;- c [Fig'l^f)]. 

The phase diagram near the BEC-QCP is obtained by min- 
imization of £'eff ^ at constant-p as a function of Vi, V2, 
and V4 (Fig. |3]). We restrict our attention to the region < 
Vi, V2, V4 <^ Jq, which is appropriate for weakly-coupled 
dimer systems with repulsive triplon-triplon interactions. Fi, 
F2, F3, and F4 are positive in this regime. Only F5 changes its 
sign on the dashed F5 = line shown in Fig. [3] The GL expan- 
sion ^ is justified because the transition is always continuous 
for the region of interest (IF5I ^ Fy, 1 < y < 4). 

The phase diagram of Fig. [3] includes eleven phases that 
are schematically represented in Fig.[3jg). (The precise wave 
functions are summarized in Table [l| in Methods.) The mag- 
netic orderings are better described in terms of S = 1/2 
pseudospin variables S/ obtained from a Matsubara-Matsuda 
transformation:^^ s^. 1/2- tj^t.^ and sf (sT) ^/i(^Ji)- 
The EEC state for H > Hd is well approximated by the mean 
field state |^) ^ n^/[cos(^,/2)^J |0), + smj2)e^^^tl |0),], 
which is a direct product of spin coherent states in terms of 
the s variables. The relations between spin and pseudospin 
expectation values are 

{St^ = (-ir^{st), (5y =^(1-2(4)), (6) 

where yu = 1 , 2 labels each spin of the dimer /. 



SINGLE-Q STATE AND ITS INSTABILITIES 

The simplest ordered phase is the single-Q state (No. 1). 
The xy pseudospin components spiral around a uniform z- 
component (cone structure). Our analysis of ^ shows that 
this phase is stable in a region where Vi, V2, and V4 are 
smaller or comparable to a characteristic energy scale 2t_i_ = 
(8/3)(/2 - -^3) ~ 3.52(2) K, where t_i_ is the intra-layer triplon 
hopping amplitude. The occupation of a single mode requires 
Fi < min(F2,F3,F4). For instance, the state No. 2, in which 
two modes ±Q„ are evenly occupied, has lower energy for 
Fi > F2. The corresponding spin configuration has a fan-like 
coplanar spatial modulation. Likewise, the condition Fi = F4 
induces an instability to the state No. 3, in which Q„ and -Q^ 
(n 4^ m) are evenly occupied. This state also exhibits the fan- 
type spatial modulation with wave-vector (Q« +Qm)/2, but the 
fan-plane rotates with a diff'erent wave- vector (Q„ - Qm)/2. 

The single-Q state also becomes unstable for Fi > F3. In 
this case, the energy is minimized by distributing the conden- 
sate over {Qi, Q2, Q3} or {-Qi, -Q2, -Q3}. FigureQc) shows 
the resulting triple-Q state (No. 4) that is a triangular super- 
structure of magnetic domains with a local approximate 120° 
xj-order. The phase diff'erence between diff'erent domains is 
+27r/3, where the ± sign distinguishes the two possible stag- 
gered chiralities of the spin superstructure. The local 120° or- 
der results from proximity of ±Q„ to the commensurate three- 
sublattice wave- vectors (weak inter-bilayer exchange) [Eq. ^ 
or Fig. [TJd)]. The lattice spacing a* of the spin superstruc- 
ture is controlled by the small shift from the commensurate 
wave- vectors: a = V3|l - 3a/(2n)\-^ ^ 23(2). Around each 



FIG. 4: Spin configurations of some ordered states: (a) No. 8 (0 = ;r, /c = +1), (b) No. 8 (0 = 0, /c = +1), (c) No. 4, (d) No. 7 (0 = n), 
and (e) No. 7 (0 = 0). To emphasize the long wave-length modulation we define the sublattice wave function if/xi^) = exp[/(AQ„ • r/ + 
Xxn)\ + exp[-/(AQ;^ • Px^^) = l^'xWI^ X = A, 5, and C, where AQ„ = [1 - 27T/(3a)]Qn is a shift of Q;^ from the closest 

wave-vector of the three- sublattice commensurate order midxxn phase factors ensuring if/xi^d = (t/> for every site / on each sublattice X. 
Cn and Cn are given in Table[l|in Methods. We show p^(r) in (a) and (c) [p^ ~ ~ Pc in (c); the inset of (c) shows Argj/r^(r)], Argif/^(r) in 
(b) (p^(r) in the inset), and Pinax(^) = niax[p^(r),p^(r),p^(r)] in (d) and (e). The large and small thick arrows in (b) indicate vortices on the 
A sublattice with vorticity -2 and +1 respectively. The double anti-vortices are located at the centers of the large dark spots in the inset, while 
the single vortices are at the center of the small dark spots. The inset of (a) illustrates the same ordered state for a small value of |Q„|. 



domain there are vortices and antivortices that form a honey- 
comb lattice [inset of Fig.|4jc)]. We will discuss more exam- 
ples of vortex crystals in the next section. 

QUANTUM VORTEX CRYSTALS 

The states so far discussed occupy less than four modes, 
meaning that the term has no contribution to the energy 
density. Our phase diagram suggests that states involving T5 
require amplitudes Vi, V2, V4 that are comparable to the ki- 
netic energy. Among such states, we next discuss the state 
(No. 8) in which the six modes are equally occupied and the 
phase parameters satisfy OJ^ = -Ikhtt/S + const, to minimize 
the energy for Fs > (see Methods |B]). Here k = ±1 is a chiral 
order parameter. This phase becomes unstable towards state 



No. 7 at the F5 = line [Figs.jSfe) andjSff)]. The new state 
is also characterized by equally occupied six modes but the 
phase variables are the same because F5 < 0: Oi = O2 = O3. 
A close examination of the wave functions of states 7, 8, and 
10 reveals an additional U(l)-invariant parameter, 

n 

that distinguishes between the spin structures shown in 
Figs. |4ja) and |4jb). The energy dependence on O is a con- 
sequence of a three-body process of the form 

~2 /]~[pQ„p_Q„cosO. (8) 
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The sign of the pref actor determines the value of 0. We will 
then discuss the two cases = 0, tt. By properly choosing r* 
(not restricted to lattice sites), the spin configuration can be 
rewritten as 



(4U 




COS 





Q„-(r, -r )-\-Sn,3- 



(9) 



Figure |4ja) shows that this 6-Q state for = tt is a crys- 
tallization of a magnetic vortices in a triangular superlat- 
tice. The crystallized vortices are not topological defects 
but thermodynami c stabl e states similar to the crystalliza- 
tions of skyrmions!^'^'^'^ Z2 vortices, and solitons,^"^ or 
the Abrikosov lattices in the type-II superconductors.^^ In 
the present case, the vortex texture appears in each of the 
three underlying sublattices (A, and C) because of proxim- 
ity of ±Q„ to the commensurate three- sublattice order wave- 
vectors. The triplon density is or nearly at the vortex cores 
and on the boundaries between adjacent vortices. The vor- 
tex crystal is chiral because all the vortices have the same 
value of K. The lattice spacing of each vortex sublattice is 
a" = V3|l - 3a/(27T)\~^. The three sublattice structure disap- 
pears for |Q„| leading to a simple triangular vortex crystal 
[see inset of Fig.Qa)]. 

Figure [4jb) illustrates the other vortex crystal obtained for 
= 0. The triangular vortex lattice has three sublattices. Two 
sublattices are occupied by single vortices (a: = ±1), while the 
third sublattice is occupied by double-antivortices (a: = =f2). 
States 5, 6, 9, 10, and 11 correspond to other vortex crystal 
structures. 

The above mentioned phase No. 7 is a 6-Q state with a 
coplanar fan-like modulation [Figs.|4jd) and|4je)]. The fan- 
plane is determined by the angle 0/2, where O = Oi = O2 = 
O3 . The spin configuration is obtained by replacing OJ^ with 
O in Eq. ^ and it is slightly diff'erent for = 0, tt. Both 
cases correspond to a triangular superlattice of magnetic do- 
mains similar to state No. 4, but the present states exhibit lo- 
cal coplanar order in each domain. The local order for = tt 
comprises three pseudospins, one of which is almost polarized 
and the other two are canted in opposite directions [Fig.|4jd)]. 
In contrast, the local order for = comprises two pseu- 
dospins canted in the same direction while the third one is 
canted in the opposite direction [Fig.|4je)]. 



EFFECTS OF ANISOTROPY 

The U(l) symmetry of our Hamiltonian for H || c is always 
broken down to the six-fold symmetry of the triangular lattice 
when exchange anisotropy induced by spin-orbit or dipole- 
dipole interactions is included. In the long wave-length limit, 
the contribution of this exchange anisotropy is 

'^-s- J](4t!q,_.2(„-1).,73^H.c) 

n 

-Z^^^^-sk-?^). (10) 



While the exchange anisotropy is small compared to the ex- 
change interactions in Ba3Mn208, this term becomes the dom- 
inant contribution very close to the QCP because it is linear in 
p, while the interaction terms included in Eq. ^ are quadratic. 
Interestingly enough, among the competing phases that appear 
in the phase diagram of Fig. |3] this term stabilizes the double- 
Q fan state (No. 2) and the vortex crystals (6-Q condensate 
No. 8 with = or = tt). *Hanis favors one of the two chiral- 
ities (/c = -1) of the the vortex crystal with net chirality. Then, 
the main role of the interaction terms ^ right above Hd is to 
split the degeneracy between these two diff'erent spin order- 
ings. The resulting phase boundary (double line in Fig.jSj is 
obtained from the condition 2(ri -h r2 - Fs - -h Fs =0. 

A diff'erent source of anisotropy appears when H is tilted 
away from the six-fold symmetry c axis because of a small 
single-ion anisotropy D ^ -0.28 K al ong the c-axis that 
breaks the six-fold symmetry into and leads to the fol- 

lowing long wave-length contribution: 



This term alone favors the coplanar states No. 2 or No. 7. The 
boundary between these phases, Fi -h F2 - F3 - F4 - F5 = 0, 
is again determined by the triplon-triplon interaction (|4]). The 
maximum amplitude of ^r^^ (H _L c) is comparable to the 
amplitude of "Hanis and the competition between these terms is 
controlled by the angle between H and the c axis. 

Thermodynamic measurements on Ba3Mn208 show two 
diff'erent phases in the vicinity of H = Hd even for H || 
A narrow Phase II that appears right above Hd and a broad 
Phase I, which is a single-Q spiral, extends over a much big- 
ger window of magnetic fields. The fact that the double line 
shown in Fig. [3] crosses the single-Q phase indicates that both 
possible scenarios could explain Phase II when the angle be- 
tween H and the c axis is small: it is either the double-Q 
fan state or a vortex crystal. This prediction can be verified 
by performing nuclear magnetic resonance (NMR) or neutron 
scattering experiments right above H = Hd- 

It is important to emphasize that our conclusions extend be- 
yond the particular case of Ba3Mn208. In general, six-fold 
symmetric frustrated magnets that exhibit six degenerate min- 
ima in the single-triplon dispersion (for quantum dimer sys- 
tems) or in the single-magnon dispersion above the satura- 
tion field (for usual magnets) should either exhibit a double-Q 
fan state or a vortex crystal close enough to the critical field. 
Therefore, the field induced magnetic vortex crystals should 
be rather ubiquitous in these materials. 



METHODS A: Computation of rq(k, k ) 

The Bethe-Salpeter equation for Fq(k, k') corresponding to 
a zero total frequency ladder diagram is 



Fq (k, W) = U-\- V, 



d3g^rq,(k,k0 (tZ+Vq-) 
StT^ a>k+q' + OJk'-q' 



. (Al) 
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Below we simplify our notation as Fq = Fq (k, k'). We will 
show that in lattice models ( |A1| ) can be reduced to coupled 
linear equations by using the following ansatz: 

rq= r ^r,. + J]A,y(r,).">-^''. (A2) 

Here denotes a displacement vector to a site within a range 
of the interaction potetntial, z is the total number of such vec- 
tors, and Aj^ = Aj^(k, W) are the coefficients to be determined. 
Since / d^q' Vq. = 0, we find 

Then, we rewrite ( |A1| ) as 

fd- 



8;r3 Wk+q' + OJk'-q' 

and introduce the following notation: 



(A3) 



(A4) 



-iq' 



By substituting ( |A2| ) into ( |A3| ), we obtain 

J]y(r,)(T'[)*A, + (To + c/-')<r> = i, 

'7=1 



Wk+q' + tt»k'-q' 

(A5) 



(A6) 



We can now saferly take the U ^ oo limit. The substitution 
of ( |A2| ) into ( |A4| ) leads to 

z 

^(T^'V(r,) + v)^v + T'[{r)= 1, l<;7<z. (A7) 



These equations can be made compact in a matrix form: 



Bn 



... Bu t\ 



"zz ' 1 



'0 J 







1 


A, 




1 









(A8) 



where Bj^y = T^Viry) + Sj^^y and C;^ = V(rj^){T^^y. Since the 
matrix elements in the left hand side are q-independent, Aj^ 
(1 < T] < z) and (F) are also q-independent as is assumed in 
the ansatz ( |A2|). While the lattice symmetry may be used to 
simplify Eg. ( |A8| ), this is a generic prescription for solving the 
integral equation ( |A1| ) for lattice models. 



METHODS B: Determination of 

The vortex lattice (No. 8) and the state No. 10 require spe- 
cial relations for 0„ (1 < ^ < 3). The relative relative phases 



are determined by the condensate densities p^Q (1 < ^ < 3) 
through the term F5 in Eq. (|4]) as is outlined below. We assume 
P+Q ^ 0. because otherwise a condition for 0„ follows in 
a trivial way. It is also easy to verify that the case with F5 < 
leads to Oi = O2 = O3. Thus we consider F5 > below. 

We can rewrite the term F5 by fixing the global U(l) phase 
at O3 = 0: 



^eff,5 = 2F5 yPQ^ P_q^Pq,P_q, 

X [cos (Oi - O2) + Ki cos Oi -h K2 cos O2] , (Bl) 



where 



PQ2P-Q2 



(B2) 



Differentiation of Seff.s with respect to <!>„ (n = 1, 2) leads to 
either 



Ao — A, Ao — A, A, — A, Ao — A^ 

cosOl = ^, COSO2 = ^ ^ ^ ^ 



or 



sinOi = sin02 = 0. 



(B3) 



(B4) 



In the first case ( |B3| ), £^eff,5 = -^5 TjhPq P_q follows, and a 
simple optimization of 

E,«=-nJ] (pQ„ + P-Q„) {pQ.. + P-Q^ 

n n 

+ {T2 - Ts) ^Pq„P_q„ + J] (pq„Pq,„ +P-q„P-q,„) 

n n<m 

+ ipQ,.P-Q,. + P-Q.PQ,,. ) (^^^ 

leads to the vortex crystal solutions. 

The other case ( |B4| ) corresponds to the state No. 10, where 
two of 0„ (1 < n < 3) are equal and the third one diff'ers 
from them by tt. Numerical analysis of Eq. ^ suggests a 
certain symmetry among p^^ for this state: p^^ = p_^^ = 

P"P^ pQ, = P_Q3 = q''p^ and Pq^ = p_Q^ = /'p (or any 
equivalent distribution obtained by a symmetry operation) as 
is given in Table |l| The eff'ective energy density is 

^eff = -2yu (2x z) + 4 (Fi -h F4) 

- 2(Fi - F2 - F3 + F4 - F5)/ + (Fi + F2)z' 

-4F5jz + 4(F3-hF4)zx, (B6) 



where x = (q'' -h /')p/2, y = Ajq''r''p, and z = p"p- By 
diff'erentiating this expression with respect to these variables 
we obtain 



-Fi + F2 + F3 - F4 + F5 



(B7) 
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TABLE I: Coefficients of the multi-Q states (t/> = exp(/Q„ • r^) + Cn exp(-/Q„ • r^)]. Equivalent states under symmetry are omitted. The 
indices of states are given in Fig.|3] For the state No. 5, /? = (Fi + Fa - 2F4)/(3Fi + F3 - 4F4) and ^ = (1 - p)/2 are required. For the state No. 9, 
/ = c(Fi+F2-F3-F4-F5)(Fi-F2+F3-F4+F5),^' = c[F2-(F3-F4)2+F2F4-Fi(F2+F4-F5)-F4F5], and / = c(Fi-F3)(Fi-F2+2F3-2F4+F5) 
are required where c is determined by + 2q' +2r' - 1. For the state No. 10, conditions for p" , q'\ and r" are obtained by solving Eqs. ( |B7| ) 
and (1^8) along with p" + q" + r" = 1/2. For the state No. 11, p'" = (Fi + F2 - F3 - F4)/(4Fi + 2F2 - 4F3 - 2F4) and q''' = 1/2 - p''' are 
required. The phase variables of states Nos. 7, 8, and 10 must satisfy = -30 + 20 + 20' = or ;7r (mod In); see Eq. ^ and the associated 
discussion. The phase variables are otherwise arbitrary. 



State 








^2 






\ 


VP 





VJ 







VJ 




VJ 


2 


-\foJ2 
















3 


















4 


4pI^ 















5 


yjpp 















6 


Vp74 














7 














8 


^ 


^^^'■^ 










9 











^/q^e'^' 




10 














11 









^|q^e''^" 








and 



ri+r4 
r3+r4 



(ri+r2)(-ri+r2+r3-r4+r5) 



-2r5 



(B8) 



and can be obtained from Eqs. ( |B7| ) and ( |B8| ) along 
with the condition p'' + q'' + = 1/2. 
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